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We study the gravitational collapse of a star with barotropic equation of state p = wp in the 
context of f(lZ) theories of gravity. Utilizing the metric formalism, we rewrite the field equations 
as those of Brans-Dicke theory with vanishing coupling parameter. By choosing the functionality of 
Ricci scalar as f(TZ) = alZ m , we show that for an appropriate initial value of the energy density, if 
a and m satisfy certain conditions, the resulting singularity would be naked, violating the cosmic 
censorship conjecture. These conditions are the ratio of the mass function to the area radius of 
the collapsing ball, negativity of the effective pressure, and the time behavior of the Kretschmann 
scalar. Also, as long as parameter a obeys certain conditions, the satisfaction of the weak energy 
condition is guaranteed by the collapsing configuration. 
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Einstein's General theory of relativity is the classical theory of one of the four fundamental forces, gravity, which 
is the weakest but most dominant force of nature governing phenomena at large scales, and is described by a math- 
ematically well-founded and elegant structure i.e., differential geometry of curved spacctime. The Einstein's held 
equations, a system of non-linear partial differential equations, relate the geometric property of spacctime to the 
four-momentum (energy density and linear momentum) of matter fields leading to precise predictions that have re- 
ceived considerable experimental confirmations with high accuracy such as solar system tests (see [l| and references 
therein). One of the most engrossing but open debates in general relativity is that of the final fate of gravitational 
collapse with possibility of the existence of spacetime singularities, the ultra-strong gravity regions where the densities 
and spacetime curvatures blow up, leading to a spacetime which is geodesically incomplete and the structure of 
any classical theory of fields is vanquished. A star with a mass many times than that of the Sun would undergo a 
continual gravitational collapse due to its self-gravity without achieving an equilibrium state in contrast to a neutron 
star or a white dwarf. Then, according to singularity theorems established by Hawking and Penrose Q a singularity 
is reached as the collapse endstate. Such a singularity may be a black hole hidden from external observers by an 
event horizon or visible to the outside Universe (naked singularity). In the latter collapse procedure, the information 
on super-dense regions can be transported via suitable non-spacelike trajectories to a distant observer. Although the 
occurrence of a spacetime singularity as the final outcome of a collapse scenario is proved by the singularity theorems, 
they do not specify the nature of such a singularity. The cosmic censorship conjecture first articulated by Penrose 
states that a black hole is always formed in complete gravitational collapse of reasonable matter fields, or a physically 
reasonable spacetime contains no naked singularities. However, up to now many exact solutions of Einstein's field 
equations describing singularities, not hidden behind an event horizon of spacetime, are known. A remarkable study 
is the one by Shapiro and Teukolsky, who showed numerically that gravitational collapse of a spheroidal dust may 
end in a naked singularity [|| . Also many exact solutions of Einstein's field equations with a variety of field-sources 
admitting naked singularities have been surveyed. The examples studied so far include gravitational collapse of a 
pressure-less matter j(| , radiation Q , perfect fluids , imperfect fluids @ and null strange quark fluids [T3| • Beside 
general theory of relativity, there exist alternative theories of gravity explaining gravitational phenomena. Such theo- 
ries have been studied for a long time [ll| . From the theoretical standpoint there has been many attempts to correct 
the Einstcin-Hilbert action in order to renormalize general relativity to build a quantum theory of gravity or at least 
some effective action (including the low-energy limit of string theories), or to quantize the scalar fields in curved 
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spacetimes [121 ] . From the observational point of view, the discovery of current acceleration of the Universe using 



CMB la supernova [13j suggests that such acceleration may be explained within the framework of general relativity by 
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assuming that 76% of energy content of the Universe is filled with a mysterious form of dark energy with equation of 
state p ~ — p (where p and p are the energy density and pressure of the cosmic fluid, respectively). Another possibility 
is to include a cosmological constant A of a very small magnitude in Einstein's field equation, but encounters such 
difficulties as the well-known cosmological constant problem and the coincidence problem. Another alteration is f(7V) 
theories of gravity [14| where the Ricci scalar in the Einstein-Hilbert Lagrangian is replaced by a general function of 
it, providing alternative gravitational models for dark energy since the explanation of the cosmic acceleration comes 
back to the fact that we do not understand gravity at large scales. Such theories can describe the transition from 
deceleration to acceleration in the evolution of the Universe fl5j . Moreover, the coincidence problem may be solved 
simply in such theories by the Universe expansion. Also some models of modified theories of gravity are predicted 
by string/M-theory considerations 0]. Recently, it has been shown that the accelerated expansion of the Universe 
may be the result of a modification to the Einstein-Hilbert action in the context of higher order gravity theories [l7| ■ 
Our purpose here is to consider the gravitational collapse of a star within the framework of f(TZ) theories of gravity, 
whose matter content obeys the barotropic equation of state, p = wp. We investigate the conditions under which 
the resulting singularity may be naked or not. In Section II we apply the metric formalism to the action of f(7V) 
gravity [lH and rewrite it as that of the Brans-Dicke theory with ujbd = 0. Choosing f(lZ) = alZ m [l9[ in Section 
III and fixing the corresponding potential, we find m as a function of initial energy density and a. In Section IV 
we study the behavior of the expansion parameter which is the key factor in examining the formation or otherwise 
of trapped surfaces during the dynamical evolution of the collapse scenario. In Section VI we examine the global 
features of the nakedness of the resulting singularity by investigating the behavior of the Kretschmann scalar as a 
function of time. In order to fully complete the model we utilize the Vaidya metric to match the interior spacetime 
to that of the exterior one. 



II. f(K) FIELD EQUATIONS 

We begin by the general action in modified theories of gravity given by 

A=— [ d 4 X^gf(K) + A ma tter(g»v,1p), (1) 

2k J 

where K = 8ttG, G is the gravitational constant, g is the determinant of the metric, TZ represents the Ricci scalar and 
ip collectively denotes the matter fields. Introducing an auxiliary field one can write the dynamically equivalent 
action as (See [2(| and references therein) 

A= A. f d^XsJ—g [/ + -¥)] + Amatterig^^), (2) 



2k _ 

where variation with respect to "J leads to the following equation as 

/"(*)(ft-¥) = 0. (3) 
If /"(*) 7^ 0, one can then recover action (JT]) by setting = 1Z. Redefining the field ^ by (p = f'{^) and setting 

V{4>) = *(M - / (4) 

action ([2]) will take the following form 

-A = 77- / d A xyf=g[4>1l - V{4>)] + A ma tter(g^,^), (5) 



2k _ 

which corresponds to the Jordan frame representation of the action of Brans-Dicke theory with Brans-Dicke parame- 
ter ujbd = 0. Brans-Dicke theory with ujbd = is sometimes called massive dilaton gravity [2l| which was originally 
suggested in [22j in order to generate a Yukawa term in the Newtonian limit. Extremizing the action yields the 
following field equations as (we set k = 8irG — 1 in the rest of this paper) 



and 



G>„=T(f), (6) 



+ 2U(0) - =T m , (7) 
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where T m stands for the trace of T™ and the subscript "to" refers to the matter fields (fields other than cj>) and we 
have defined the effective stress-energy tensor as 

Tjf) = i(T-+T^), (8) 

with 

T*, = (V^V^ - g^U<t>) - ^V(<fi), (9) 

and 

T^ m = diag (p m ,p m ,p m ,Pm) , (fO) 
being the stress-energy tensors of the scalar field and a perfect fluid, respectively. 



III. GRAVITATIONAL COLLAPSE OF A HOMOGENEOUS CLOUD WITH f(TZ) = a1Z 



Let us now build and investigate a homogeneous class of collapsing models in f(lZ) gravity with to ^ 0, where the 
trapping of light is avoided till the formation of singularity, allowing the singularity to be visible to outside observers. 
In order to achieve our purpose we examine a spherically symmetric homogeneous scalar field, 4> = $(r) originating 
from geometry. Since the interior spacetime is a dynamical one, we parameterize its line element as follows 

ds 2 = -dr 2 + a 2 {T){dr 2 + r 2 dn 2 ), (11) 

where r is the proper time of a free falling observer whose geodesic trajectories are distinguished by the comoving 
radial coordinate r and dfl 2 is the standard line element on the unit 2-sphere. It is worth mentioning that here we 
assume that starting from the homogeneous initial data, the collapsing configuration remains homogeneous till the 
singularity is formed. But as the collapse proceeds there may be some inhomogeneities occurring throughout the 
collapse scenario, the existence of which can be investigated by perturbation theory, that is, imposing inhomogeneous 
perturbations on the energy density, scale factor and BD scalar field and then see whether the terms rising from 
inhomogeneity are dominant in the formation of the singularity or not. Here we do not deal with such an issue but 
for more details the reader may consult (42j and references there in. Since the presence of matter acting as a "seed" 
field prompts the collapse of the BD scalar field, we have considered perfect fluid models obeying barotropic equation 
of state as 



WPn 



(12) 



Using the conservation equation for the matter (V Q T^ = 0) together with the use of above equation, one gets the 
following relations between p m , p m and the scale factor as 



Pm = p n , n a- 3{1+w) ; Pm = wp, n a^ 1+w \ 



(13) 



where p Qrn — p m (a = 1), is the initial value of energy density of matter on the collapsing volume. Making use of 
equation ([8]) and equation (|11[) one finds the following equations for the effective stress-energy tensor as 



P^-- TT[ r S) - l -{Prn+P,)= 1 - 



A 

a 



v(4>) 



(14) 



and 



p (cff) = T r ^ = T s (° ff) = T v > = I (p m + p,) = i 



2-c 



V{4>) 



(15) 



with all other off-diagonal terms being zero and the radial and tangential profiles of pressure are equal due to the 
homogeneity and isotropy. Substituting the line element ()lll) into Einstein's equation one gets the interior solution as 



M' 
R 2 R' 



P(cff) 



M 

R 2 R' 



(16) 
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R -T' 



(17) 



where A4(r, r) rises as a free function from the integration of Einstein's field equation which can be interpreted 
physically as the total mass within the collapsing cloud at a coordinate radius r with M. > 0, and R(r,r) = ra(r) is 
the physical area radius for the volume labeled by the comoving coordinate r. From equation (|14j) and first part of 
equation (fH)]) , one can solve for the mass function 



R 3 



(18) 



Using the above equation together with equation (|17p we arrive at a relation between a and the effective energy 
density as follows 



(P 4 +Pm)- 



(19) 



Since we are concerned with a continual collapsing scenario, the time derivative of the scale factor should be negative 
(a < 0) implying that the physical area radius of the collapsing volume for constant value of r decreases monotonically. 
The singularity arising as the final state of collapse at r = r s is given by a(r s ) = 0. On the other hand when the 
scale factor and physical area radius of all the collapsing shells vanish, the collapsing cloud has reached a singularity. 
A point at which the energy density and pressure blows up, the Kretschmann scalar JC = R abcd R a bcd diverges and 
the normal differentiability and manifold structures break down. In order to solve the field equations we proceed by 
substituting for /O (off) and P (off) from equation (|16p into equations (|14[) and (| 15[) and rewrite them as follows 



1 



(20) 



- +2- 



1 



wp TO a" 3(1+Iu) + am$ + 2am-$ - 



(21) 



where for later convenience we have rescaled the scalar field as <j> — am$ (a and m are real constants), and by the 
virtue of equation (j4|) the associated potential to f(TZ) = a7Z m can be fixed as 



V(&) = a(m — l)$m-i ; m^l. 



(22) 



Since the scalar field must diverge at the singularity we examine its behavior by taking the following ansatz for the 
scalar field 



$(r) = a 5 (r), 



(23) 



where 6 is a constant whose sign decides the divergence of the scalar field. Substituting the first and second time 
derivatives of the scale factor from equations ([20]) and ([21]) into equation {7J together with the use of equation (|22]) 
one finds 



-[5+3(1 



i +w}] f2p om (8 + 3w-l) \ r_y f 4+j(2m-l) 



3am5(2 + S) 



2m 



6mS + 3m<5 2 



Matching the powers of scale factor in equation (j2~4")l we arrive at the following expression for 5 



6 = 



3(1 + w)(l - m) 



(24) 



(25) 



whence by substituting the above equation into the pair of square brackets in equation (|24[) one gets the following 
equation to be satisfied by m 



a{3(l + to) + m [-13 - 9w + m(8 + 6w)} | + 2p om [Am - 3(1 + to)] 



0. 



(26) 
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Solving the above equation, we find m as a function of a and p m for w as 



for w = 0, 



13a 



± 



m± = 



73a 2 - 16ap m + 64p 2 



16a 



(27) 



for w 



m± 



5a + Ap Qm ± 13a 2 - Wap om + 16p 2 



Go 



(28) 



for w 



■|, and 



-7a + 8p m ± 



33a 2 - 80ap om + 64/r^ 



8o 



(29) 



for to = i, where a^O. 



m± = 



4a - 2p om ±^2 3a 2 + 2ap QWi + 2p 2 n 



5a 



(30) 



IV. TIME BEHAVIOR OF THE SCALE FACTOR AND SINGULAR EPOCH 

One would like to study time behavior of the scale factor as the collapse evolves, considering matter fields. If at 
time t — t* (or cquivalcntly for some a = a*) the collapse begins, then by integrating equation (|20|) together with 
the use of equations (|2"2"|) and (|2"3"j) near the singularity with respect to time one gets the time behavior of the scale 
factor as 



a(r) 



i,,.,,, . 1 2p + aim — 1) 



5 + 3(l + m) 



and the corresponding singular epoch as 



T. =2\ 



l2p (>m + a(m - 1) a *i(S+3(i+w)) 
6am(l + 8) S + 3(l+w) 



(31) 



(32) 



where the time t s corresponds to a vanishing scale factor. Thus the collapse reaches the singularity in a finite proper 
time. This result for the scale factor completes the interior solution within the collapsing cloud, providing us with 
the required construction. 



V. THE CONDITIONS 



We are now in a position to investigate the nature of singularity as the endstatc of a collapsing scenario. The 
singularity is called locally naked if it is only visible to an observer being in the neighborhood of it (such a singularity 
is necessarily covered by a spacetime event horizon) and is called globally naked if there exists a family of future 
directed non-spacelike geodesies reaching to the outside observers in the spacetime and terminating at the past in the 
singularity. To investigate the nature of spacetime singularity arising from the collapse procedure, we examine here 
whether such a singularity could be naked, or necessarily covered within a spacetime event horizon and if so under 
what conditions. Since a can be regarded as a free parameter, we seek for the appropriate values of it satisfying the 
following conditions 
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• The ratio Ai/R stays less than unity till the singular epoch which means that the singularity has been formed 
earlier than the formation of the apparent horizon or the formation of trapped surfaces have been failed till the 
singular epoch. 



• S < during the gravitational collapse scenario accompanied by divergence of the scalar field. 



• For physical reason, weak energy condition, stated as p (off) > and p (off) + P (off) > must be satisfied during 
the dynamical evolution of the system. 



• The effective energy density and effective pressure blow up in the vicinity of the singularity, the latter being 
negative during the evolution of the collapse process, since the absence of trapped surfaces is accompanied by 
a negative pressure. In fact the negativeness of the effective pressure implies that Ai < 0, that is, the mass 
contained in the collapsing volume with comoving coordinate r keeps decreasing leading to an outward energy 
flux during the gravitational collapse scenario. 



Kretschmann scalar diverges at the singular time and then converges to zero at late times. 



In order to determine whether the singularity is naked or not, one needs to investigate the formation of trapped surfaces 
during the collapse procedure. These surfaces are defined as compact two-dimensional (smooth) spacclike surfaces such 
that both families of ingoing and outgoing null geodesies orthogonal to them necessarily converge or the expansion 
parameter 9 of the outgoing future-directed null geodesies is everywhere negative [23] ], Consider a congruence of 
outgoing radial null geodesies having the tangent vector (V T , V r , 0, 0), where V T — dr/dk and V r = dr/dk and k is 
an affine parameter along the geodesies. For the spacctimc metric (II ID . the geodesic expansion parameter which is 
defined as the covariant divergence of the vector field V v is given by [24[ 



9 = W V V V = - 
r 



V 



(33) 



If the null geodesies terminate at the singularity in the past with a definite tangent, then at the singularity we have 
> 0. If such family of curves do not exist and the event horizon forms earlier than the singularity, a black hole is 
formed. Utilizing equation p3[) . we now attempt to study the formation of trapped surfaces during the dynamical 
evolution of the gravitational collapse procedure. We show that physically, the formation of a black hole or a naked 
singularity as the final state for the dynamical evolution is governed by the rate of collapse and the presence of scalar 
field. It is seen that for a specified range of variation of a, the cosmic censorship conjecture may be violated for all 
cases of matter considered below. In the following subsections we consider first the four conditions mentioned in the 
beginning of this section and postpone the last one to the next section. We begin by calculating the ratio Ai/R in 
the general case which is considered for the four cases of matter, w = {0, — |, — |, |} corresponding to dust, cosmic 
strings, domain walls and radiation, respectively. By the virtue of equation (|17|) we have 



M 
R 



2p n 



a(m — I) 



6am(l + 5) 



-(S+3w + l) 



(34) 



The weak energy condition which states that the energy density as measured by any local observer must be non- 
negative can be written for any timelike vector as follows 



T^WV > 0, 

whereby one gets the following conditions for the effective energy density (p (off) > 0) 



P(off) 



2p m + a(m — 1) 



2am(l + 5) 

and the sum of effective energy density and pressure (/O (off) + p (off) > 0) as 

2p om + a(m — 1) 
2am(l + 5) 



(<5+3(™+i)) > q 



P(cff) +2V) = I l + t» + 



-(5+3(«>+l)) 



> 0. 



(35) 



(36) 



(37) 
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Finally for the rate of change of mass function with respect to time one has 



M = r 3 (S + 3w) 



2p m + a(m — 1) 
6am(l + S) 



a -(S+3 W +l) | h 



(38) 



where the minus sign has been absorbed into a. At the initial epoch where a(r*) = 1 there should not be any trapping 
of light, then by assuming that r = is the boundary of the collapsing volume one may easily see that for a suitable 
initial value of energy density the ratio Ai / R is less than unity at the initial time. This fact is in accordance with 
the regularity conditions stating that if gravitationally collapsing massive stars are to be modeled, then the energy 
density, pressure, and other physical quantities must be finite and regular at the initial spacelike hyper-surface from 
which the collapse commences. For the case of homogeneous-density collapse the resulting singularity coincides with 
the curves R(t s , 0) = or R(t s , r ^ 0) = 0, corresponding to a central or non-central singularity, respectively. In the 
next subsections we first consider the simpler case of non-central singularity and investigate formation or otherwise 
of trapped surfaces for different values of w. 



A. Dust (w = 0) 

For this case of matter we have the following relations (we set p = 1 in the rest of this paper) 



M _ r 2 
R 1 



2 P,im+ a ( m_1 ) 



6am( 1+5) 



-(5+1) 



P(o ff) +P (eB) = (! + !) 



ZPjjm+aCm-l) 



2am(l+5) 



2am(l+<5) 



"(-5+3) 



Ptefi) 



M = r 3 6 



-(5+3) 



2p nm +a("»-l) 



2am(l+5) 



2 P nm 



6am(l+<5) 



-(5+3) 



,-(5+1) 



(39) 



If the interval for a is —3.9 < a < —0.1 and by picking out m_ in equation (|27j) . the corresponding parameter 
<5_ is always negative, its absolute value varies between 1 <| 5- \< 3 and its minimum and maximum values are 
5- = —2.72484 and 6- = —1.00569, respectively. It is seen that for such values of a and the ratio A4/R stays 
less than unity and then the expansion parameter is always positive up to the singularity, that is the singularity is 
formed earlier than the formation of apparent horizon which is the boundary of trapped surfaces. The negativeness 
of effective pressure for the allowed values of a ensures that M. is negative as collapse proceeds which means that 
the mass contained in the collapsing volume keeps waning. Then there exists an outward energy flux which may be 
visible to outside observers (for the case of globally naked singularity) since the trapped surfaces do not form early 
enough to cover the singularity. In addition, the validity of the weak energy condition can be easily checked by the 
virtue of expressions obtained for p (off) and p (cff ) + p (off) . 



B. Cosmic Strings (w = — |) 

Cosmic strings are the result of hypothetical 1-dimensional (spatially) topological defects which may have been 
constructed during a symmetry breaking phase transition at the early Universe. The possibility of their existence 
was first considered by Tom Kibble in 1976 J2^|. A fluid of cosmic strings may have an effective equation of state, 
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hp , so one has the following relations for this type of matter fluid as 



M 



2 P nm + Q ( m - 1 ) 



6am(l+<5) 



p (oft) +p (off ) 



'5+21 
> 3 ; 



Pf. 



' 5-1 
, 3 



2« +a(m-l) 



2am(l+j) 



2am(l+5) 



-(5+2) 



2 P Qm +"('"- 1 ) 

2aro(l+5) 



.M = r 3 (£-l) 



-(5+2) 



ZPjjm+aCm-l) 



6am(l+<5) 



-(5+2) 



(40) 



a- 5 I d 



Choosing m + in equation (|28j) , it can be seen that the valid range of change for a is —100 < a < —0.1 upon which the 
absolute value of parameter 8 + is restricted to vary between <| <5+ |< 2 (in order to prevent the ratio M./R to be 
infinity, 8 + = —1 has to be excluded) and its minimum and maximum values are 8+ = —1.86224 and 8+ = —0.61558, 
respectively. It should be noted that one can choose the lower limit of a to be much less than —100, but such a 
choice does not affect considerably the value of 6+ and its magnitude remains close to —0.6. Taking these values into 
account, the ratio M. /R stays finite till the singular epoch and causes the expansion parameter to be positive up to 
the singularity, and if no trapped surfaces exist initially then none would form until the epoch a(r s ) = which is 
consistent with the fact that there exist families of outgoing radial null geodesies emerging from the singularity. Also 
the weak energy condition is satisfied (| 6+ |< 2) and the effective pressure is negative (since 8+ < 0), consistent with 
the fact that time derivative of mass function is negative i.e., the mass incorporated in the region where the collapse 
procedure evolves keeps falling off. As a result, there exists an outward energy flux during the collapse scenario which 
may be visible to external Universe. 



C. Domain Walls (w = -§) 



Domain walls are two-dimensional objects that form when a discrete symmetry is spontaneously broken at a phase 
transition l26l . It has been noticed that there exist a link between domain-walls and cosmologies such as branc 
cosmology [27|. The effective equation of state for a fluid of domain walls may be p = — \p m and the mentioned 
conditions for this case can be written as 



M 
R 



P(efl) 



2p m +a(m— 1) 



6am(l+<5) 



,-(s-i) 



< P(e«)+P( 8ff , = (^) 



' 5-2 ' 
- 3 - 



2 P nm + a (™- 1 ) 



2am(l+5) 



-(5+1) 



2p Q _ m +a(™-l) 



M = r 3 (S-2) 



2p nm +a(m-l) 



2am(l+5) 



-(5+1) 



2am(l+5) 



2p nTO +a(m-l) 



6am(l+i5) 



,-(5+1) _ 



,.-(<5-D 



(41) 



As long as a varies in the range —100 < a < —0.1 (lower limit of a can be chosen much less than —100 but such 
a choice has no noticeable influence on the maximum value of 8), by choosing m + in equation (|29[) . the parameter 
8+ is negative, its absolute value is less than unity and its minimum and maximum values are 5+ = —0.953501 and 
8 + = —0.379572, respectively. One then may easily see that at initial epoch (a(r*) = 1), the regularity condition 
(there should be no trapped surfaces at the initial hyper-surface from which the collapse commences) is satisfied 
and the ratio of mass function to physical area radius of the collapsing volume is less than unity during the collapse 
procedure denoting that the expansion parameter being positive up to the singularity. In this case the collapse to a 
naked singularity may take place, where the trapped surfaces do not form early enough or arc avoided in the spacetimc. 
Also the mass contained in the collapsing ball reduces as the time advances due to the fact that the effective pressure 
stays negative till the singular epoch. It is obvious that for such values of a and 8+ the weak energy condition is 
satisfied during the collapse scenario. 
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D. Radiation (w — |) 

Finally, for this type of matter we have following relations for the said conditions 

-2p nm +a(m-iy 



M _ 2 
R — 1 



6am(l+5) 



P(off) +-P(off) { 1} ) 



2am(l+S) 



-(5+4) 



P(eff) 



2» m +a(m-l) 



M = r 3 (S+ 1) 



2am(l+5) 



2o +a(m-l) 



6am(l+i5) 



-(5+4) 



(42) 



Pf. 



2am(l+{) 



-(5+4) 



The suitable range of variation of a for m_ in equation (|30p is —0.66 < a < —0.01 which makes the corresponding 
values of 6- to vary in the range 2 <| 6- |< 4 with 6- = —2.01007 and 6- = —3.95037 are the maximum and 
minimum values of this parameter respectively. For such values of a and <5_ the ratio M./R stays less than unity as 
the collapse procedure ends. Then the expansion parameter remains positive up to the singularity denoting that the 
apparent horizon is failed to form. On the other hand, weak energy condition is satisfied for these values of a and (5_ 
and the effective pressure is negative till the singular epoch. It is worth noticing that if one chooses the parameter a 
and the power of Ricci scalar in equations (f27| -(|30 |) other than the ones determined in the above subsections, weak 
energy condition together with the conditions on effective pressure and the ratio M/R would be violated. The central 
singularity occurring at R = 0, r = can be naked if we have any future-directed outgoing null geodesies terminating 
in the past at the singularity. In order to examine the possibility of existence of such families let us introduce a new 
variable y = r 1 with 7 > 1 defined in such away that the ratio i?'/r 7_1 is a unique finite quantity in the limit r — > 0. 
Now consider the outgoing radial null geodesic equation which is given by 



dr 

— = o(r). 
dr 



In terms of variables y = r 7 and R the above equation reads 



dR 
dy 



1 



1 >' 



7-1 



^dr . 
R—+R' 

dr 



(43) 



(44) 



whence using equation (TT7)) we have 



dR 
dy 



R' 



(45) 



If there exist outgoing radial null geodesies in the past at the central singularity which occurs at t = r s , then along 
such geodesies we have R — > as r — > 0, or in terms of the variables y and R, the point y = 0, R = is a singularity 
of the above first order differential equation. For such a congruence of geodesies dR/dy must be positive up to the 
singularity. Then as long as a, p m , and m satisfy the values which was determined previously, trapped surfaces would 
fail to form as the collapse evolves. 



VI. NAKEDNESS OF THE SINGULARITY 



The final fate of a continual gravitational collapse of a matter cloud ends in either a black hole or a naked singularity 
where in the former there exists an event horizon of spacetime developing earlier than the formation of the singularity 
to cover it. Thus the regions of extreme curvatures and densities are concealed from the outside observers. The 
event horizon or surface of a black hole is defined as the boundary of the spacetime that is causally connected to 
future null infinity (28| . or in other word the boundary between events from which light rays emitted inside this 
boundary surface can not escape to future infinity while those emitted outside in a suitable direction can. Both event 
and apparent horizons coincide in stationary spacetimes, however this case is not generally true in dynamical ones. 
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Although the existence of an apparent horizon predicates the existence of spacetime event horizon, the converse is 
not always true and the event horizon may veil the singularity even if apparent horizon does not emerge on a spatial 
slice. So far we have discussed the conditions under which the collapse scenario ends in a locally naked singularity, 
that is, the singularity is visible to an observer being in the neighborhood of it. In this case, the trajectories coming 
out of the singularity do not actually come out to a distant observer but fall back into the singularity again at a 
later time without going out of the boundary of the star. Thus the locally naked singularity could still be covered 
by the event horizon and only strong version of cosmic censorship conjecture is violated but the week form of it is 
intact 29] . However if the event horizon is delayed to form or the singularity forms early enough before the formation 
of event horizon, then it would be visible to external observers and thus a globally naked singularity would born 
as the endstate of collapse rather than a black hole. Therefore in such a situation curvature invariants namely the 
Kretschmann scalar should increase near the singularity, diverge at the singular epoch and then converge to zero at 
late times [30| . In previous section we showed that for suitable values of a, trapping of light can be avoided which 
means that the apparent horizon is failed to form. But since the failure of formation of an apparent horizon does 
not necessarily bode the absence of an event horizon, we investigate the nakedness of the singularity in spherically 
symmetric collapse of a fluid by considering the behavior of Kretschmann scalar with respect to time. For the line 
clement (|11[) this quantity is given by 

K = R ahcd Rabcd = — A [a 2 a 2 + a 4 ] . (46) 

By the virtue of equation pip one can easily obtain this quantity as a function of time for w = {0, — 5, — |, ^} and 
the results are sketched in Figures 1-4. It is seen that Kretschmann scalar diverges at singular time and then tends to 
zero at late times. In order to better understand the situation one should resort to critical behavior at the black hole 
threshold. Such a behavior was discovered in gravitational collapse of a spherically symmetric massless scalar field j3l| , 



axisymmetric gravitational waves [32j, spherical system of a radiation fluid 33[ and spherical system of a perfect fluid 
obeying the equation of state p = wp J34|. Consider a type II critical collapse in which the black hole mass is scaled 
as Mbh oc\ p~p* I 7 j where p parameterizes a family of initial data sets evolving through Einstein's equations, is a 



critical value and 7 is a positive constant which is called a critical exponent [30,[35|. For a sufficiently large value of 



p the collapse procedure develops to a black hole and for a sufficiently small one it evolves to a dispersion [35J. The 
boundary between these two regimes is the black hole threshold. Now consider the limit from supercritical collapse 
(p > p*) to a critical collapse, i.e., p — > p*. In such a limit, the black hole mass tends to zero and the maximum 
value of curvature diverges just outside the event horizon. Since we have arbitrarily strong curvature outside the 
event horizon by fine-tuning, the black hole threshold can be regarded as a globally naked singularity [Hj|. Let us 
now consider the geometry of the exterior spacetime. In order to complete the model we need to match the interior 
spacetime of the dynamical collapse to a suitable exterior geometry. The Schwarzschild solution is a useful model to 
describe the spacetime outside stars but the spacetime outside such a star may be filled with radiated energy from 
the star in the form of electromagnetic radiation. The Schwarzschild solution does not describe this properly as it 
deals with an empty spacetime given by T a & = 0. The spacetime outside a spherically symmetric star surrounded by 



radiation emitted from the star is described by the Vaidya metric [36| which can be given in the following form as 

2M(r v ,v)~ 



dSout = - 



1 - 



dv z ~2dvdr v +r z v dn\ (47) 



where v is the retarded null coordinate, r v and M(r Vl v) are the Vaidya radius and Vaidya mass, respectively. In what 
follows, we use the Isreal-Darmois junction conditions to match the interior spacetime to a Vaidya exterior geometry 
at the boundary hyper-surface S given by r — r b . The spacetime metric just inside S is given by 

ds 2 n = -dr 2 + a 2 (r) [dr 2 + r 2 b dtt 2 } , (48) 

whereby matching the area radius at the boundary one gets 

R(r b ,T)=r v (v). (49) 

One then gets the interior and exterior metrics on the hyper-surface E as follows 

dsl in = -dr 2 + a 2 (r)r 2 dn 2 1 (50) 



d« 2 

as T,out 



1 _ 2M(r v ,v) | 2 dr v 
r v dv 



dv 2 +r 2 v dn 2 , (51) 
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where matching the first fundamental form gives 



dv 



1 



2M(r„«) 



dv 



(r v )s = r b a{r). 



(52) 



In order to match the second fundamental form (extrinsic curvature) for interior and exterior spacetimes we need to 
find the unit vector field normal to the hyper-surface E. We then proceed by taking into account the fact that any 
spacetime metric can be written locally in the following form as 



ds 2 



(a 2 - A/3 1 ) dr 2 - 2^dx l dr + h^dx 1 , 



(53) 



where a, f3 l , and h t j are the lapse function, shift vector, and induced metric, respectively and i, j run in {1,2,3}. 
Comparing the above equation with equations (|47[) and (|48|) together with using the following normalization condition 
for n v and n Tv 



n v n v + n rv n r ^ = 1, 
one gets the normal vector fields for the interior and exterior spacetimes as 

<„ = [0,a(r)- 1 ,0,0], 



(54) 
(55) 



1 - 



2M(r z ,,v) 



dr v 
dv 



l 



2tf(r„ii) 



2 dr^ 
dv 



l - 



2M(r v ,v) 



2^ 

dv 



(56) 



The extrinsic curvature of the hyper-surface E is defined as the Lie derivative of the metric tensor with respect to the 
normal vector n, given by the following relation as 



Kab = \Cn9ab = - [g a b, c n c + g c bn^ a + g ac rf b ] , 



2 2 

whereby the nonzero components of the extrinsic curvature read 

1 - 



oo 



r b a{T), 



K' 



r v dv 



_ 2M(r v ,v) , 2 dr v 
Tn, dv 



(57) 



(58) 



Setting [Kgg — Kg/}*] = on the hyper-surface E, together by using equations (fTT)) and ([52^1 one gets the following 
relation between the mass function and Vaidya mass on the boundary as 



M(T,r b ) = 2M(r v ,v). 



(59) 



From equation (|59)) and equation ([18]) it is seen that the behavior of Vaidya mass is decided by the allowed values 
of p m , ol and m which prompt the gravitational collapse scenario to end in the formation of a naked singularity. 
In order to get a relation describing the rate of change of Vaidya mass with respect to r v one has to match the r 
component of the extrinsic curvature on the hyper-surface E. Setting \K™ T — K™*] =0we have 



M{r v ,v) >rv = — + rla'd. 



(60) 



Now it can be seen that at the singular time, r = t s the ratio 2M(r v , v)/r v tends to zero. Thus the exterior spacetime 
at the singular epoch reads 



ds A 



-dv 2 



2dvdr„ + rldri 2 



(61) 



which describes a Minkowski spacetime in retarded null coordinates. Hence, the exterior generalized Vaidya metric 
at singular time can be smoothly extended to the Minkowski spacetime as the collapse completes [37|. The occurrence 
of a naked singularity as the final fate of a collapse scenario depends on the existence of families of non-spacelike 
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trajectories reaching faraway observers and terminating in the past at the singularity. In order to show this we begin 
by equation (f52"j) and after using equation (f5!)|) we get 



dv 



1 - r b a 



1 



M{r,r b ) 



(62) 



It is seen that imposing the null condition on the Vaidya metric leads to the same relation as the above. This means 
that null geodesies can come out from the singularity and reach distant observers before it evaporates into the free 
space. On the other hand, since for the allowed values of a, p m and m formation of trapped surfaces in spacetime is 
avoided and from another side the singularity emerges outside of the event horizon, such a congruence of trajectories 
can be detected by the outside observer. 



VII. CONCLUSION AND OUTLOOK 



One of the physical motivations for discussing naked singularities is that these objects provide a useful laboratory 
for quantum gravity, since in such ultra-strong gravity regions the length and time scales are comparable to the Planck 
length and time. In other words, quantum effects occurring in such super-dense regimes are no longer covered by the 
spacetime event horizon and the chance to observe such effects in the universe is provided. An example is quantum 
particle creation due to the formation of a naked singularity, which has been studied in the literature |38( . During the 
past twenty years, cosmic censorship conjecture has been extensively investigated in spherical models of gravitational 
collapse of physically reasonable matter. The simplest of those which has been scrutinized in detail and has been 
shown that both black holes and naked singularities form from generic initial conditions is gravitational collapse of a 
dust fluid. Dwivcdi and Joshi in (39j and Waugh and Lake in [40| showed that a naked strong curvature singularity 
can be formed in an inhomogeneous dust collapse and a self-similar one, respectively. Also some examples of naked 
singularity formation in gravitational collapse of a scalar field is given in [4l|. In this work we have studied the 
process of gravitational collapse of a star where the matter fluid obeys the barotropic equation of state p = wp, in 
the context of f(7V) theories of gravity. Making use of metric formalism we wrote the action of f(lZ) gravity as the 
Brans-Dicke one with vanishing coupling parameter. Having solved the resulting field equations by taking the ansatz 
(f2Ul) for scalar field we arrived at the expressions (|27|) - (f5U|) for the exponent of Ricci scalar as a function of a and 
initial energy density. In Section V we imposed five conditions on the effective energy density and pressure, the ratio 
Ai I R, parameter <5, time behavior of the mass function, and Krctschmann scalar, the validity of which depends on 
determining appropriate values of a. As long as these conditions are fulfilled the resulting singularity can be globally 
naked, i.e., ultra-dense regions are no longer covered by a spacetime event horizon and physical effects are allowed 
to be shared by the external Universe. It is worth mentioning that there are future finite-time singularities in the 
dark energy universe coming from modified gravity as well as in other dark energy theories. Considering f(lZ) gravity 
models that satisfy cosmological viability conditions (chameleon mechanism), it is possible to show that finite-time 
singularities emerge in several cases. Such singularities can be classified according to the values of the scale factor 
a(t), the density p. and the pressure p [Hj]. 
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FIG. 1: The behavior of Kretschmann scalar (in units of s 4 ) as a function of proper time for w = A (upper-left figure) and different 
values of aandL. a = -0.65, m_ = 2.02575 and 6- = -2.02542, Solid curve, a = -0.15, m_ = 5.97355 and 5_ = -3.33038, Dotted 
curve and a = —0.05, m_ = 6.6128 and <5_ = —3.75922, Dashed curve. For the initial energy density, scale factor, and proper time we 
have adopted the values p 0m = 1, a* = 1, and t* = 0, respectively. The corresponding singular epoches are, t b = 2.49674 for Solid curve, 
r s = 9.44081 for Dotted curve, and r s = 27.8945 for Dashed curve. 

FIG. 2: The behavior of Kretschmann scalar (in units of s~ 4 ) as a function of proper time for w = (upper-right figure) and different 
values of a and <5_. a = -3.5, m_ = 1.52406 and <5_ = -1.03175, Solid curve, a = -2.5, m_ = 1.60424 and <5_ = -1.12996, Dotted 
curve and a = —1.5, m_ = 1.8076 and <5_ = —1.34034, Dashed curve. For the initial energy density, scale factor, and proper time we 
have adopted the values p 0m = 1, a* = 1, and r* = 0, respectively. The corresponding singular epoches are, r B = 2.5083 for Solid curve, 
t s = 2.70354 for Dotted curve, and t s = 3.19309 for Dashed curve. 

FIG. 3: The behavior of Kretschmann scalar (in units of s -4 ) as a function of proper time for w = — A (lower-left figure) and different 
values of a and S+. a = -90, m+ = 1.44581 and 5+ = -0.616689, Solid curve, a = -30, m + = 1.46909 and 8+ = -0.638609, Dotted 
curve and a = —10, m+ = 1.54031 and (5+ = —0.701562, Dashed curve. For the initial energy density, scale factor, and proper time we 
have adopted the values p 0m = 1, a* = 1, and r* = 0, respectively. The corresponding singular epoches are, t b = 4.04176 for Solid curve, 
t s = 4.13327 for Dotted curve, and r s = 4.38509 for Dashed curve. 



FIG. 4: The behavior of Kretschmann scalar (in units of s 4 ) as a function of proper time for w = — | (lower-right figure) and different 
values of a and S+. a = -100, m+ = 1.61179 and S+ = -0.379572, Solid curve, a = -20, m+ = 1.68699 and <5 + = -0.407227, Dotted 
curve and a = —2, m+ = 2.55425 and <5_|_ = —0.608495, Dashed curve. For the initial energy density, scale factor, and proper time we 
have adopted the values p 0m = 1, a* = 1, and t* = 0, respectively. The corresponding singular epoches are, r s = 10.2643 for Solid curve, 
t s = 10.787 for Dotted curve, and t s = 16.808 for Dashed curve. 
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